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CHARACTERISTICS  OF  IMAGES  DERIVED 
FROM  FACTORABLE  POLYNOMIALS 


INTRODUCTION 


In  recent  years  image  processing  techniques  have 
been  developed  which  use  as  an  input  the  autocorrela¬ 
tion  function  of  an  image  or,  equivalently,  the  modulus 
square  of  the  Fourier  transform  of  the  image. 
Fierup's^  technique  is  a  good  example  of  this  approach. 
In  his  technique  an  iteration  procedure  is  used  which 
alternates  between  spatial  and  frequency  domains, 
enforcing  known  constraints  in  each  domain.  In  the 
spatial  domain  the  constraints  are  ror.-r  egativ  ity  of 
the  image  brightness  values  and  limited  spatial  extent. 
In  the  frequency  domain  the  constraint  is  that  the  pro¬ 
cessed  image  Fourier  transform  must  have  a  modulus 
equal  to  the  original  modulus.  This  processing  has  had 
a  premising  degree  of  success  on  the  limited  number  of 
images  processed  thus  far. 


In  the  absence  of  the  ren-regativity  constraint 
there  are  an  infinite  number  of  objects  which  can  pro¬ 
duce  the  same  autocorrelation  function  since  the  phase 
component  of  the  Fourier  transform  can  be  varied  at 
will  without  affecting  the  image  autocorrelation  func¬ 
tion  while  affecting  the  image  significantly,  usually 


creating  negative  image  values.  The  enforcement  of 
non-negativity  of  the  image,  however,  limits  the  number 
of  images  which  have  the  same  autocorrelation  function. 

An  important  question  is  how  many  non-negative 
images  can  produce  the  same  autocorrelation  function? 
There  seems  to  be  no  general  theory  for  two  dimensional 
images.  There  has  beer,  some  work  done  or*,  specialized 
two  dimensional  images.  Fried^,  following  the  work  of 
Bruck  and  Sodin^,  has  developed  a  method  for  developing 
specialized  two  dimensional  images  which  have  a  high 
degree  of  ambiguity  ir.  the  sense  that  a  large  number  of 
images  can  produce  the  same  autocorrelation  function 
and  satisfy  the  non-negative  constraint.  These  are 
images  derived  from  completely  factorable  polynomial 
functions.  Left  unanswered  in  Fried's  work  was  the 
question  of  whether  such  specialized  images  are  of  any 
interest.  Are  they  similar  to  real  world  images  of 
interest  and  are  the  ambiguous  images  different  from 
each  other  in  an  interesting  way?  This  report 
addresses  these  questions. 
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REVIEW  OF  THEORY 

Let  Apq  represent  the  brightness  value  of  an  image 
in  p,  q  coordinates.  We  define  a  polynomial  function 
f(u,v)  as 

f(u,v)  =  5  A  upvq.  (1) 

pq 

If  the  polynomial  has  the  characteristics  of  being  com¬ 
pletely  factorable,  then  f(u,v)  can  also  be  written  as 

f ( u  ,  v )  =  TT  (u  +  ar)(v  +  bs).  (2) 

r  ,s 

In  Fried's  and  Bruck  and  Sodin's  work  the  form  used  was 
(u  -  ar).  The  change  of  sign  here  is  for  bookkeeping 
convenience  in  computer  programs. 
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The  general  image  A  will  not  have  the  charac- 

P 

teristic  of  being  completely  factorable.  We  are  deal¬ 
ing  here  with  special  cases.  Assuming  the  image  is 
completely  factorable  then  Equations  (1)  and  (2)  can  be 
equated . 
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pq 


pq 


uPvq  ■ir  (  u  +  ar )  ( ' 
r  ,s 


bs> 


(3) 


Given  the  set  of  (a  h  )  the  image  values  A__  are  found 

r  f  5  pq 

by  multiplying  the  factors  on  the  right  hand  side  of 
Eq .  3 »  and  equating  the  coefficients  of  like  powers  of  u 
and  v.  For  example  let  the  factorable  polynomials  be 

f(u,v)  =  (u  +1  )(u  +  3)(v  +  2 ) ( v  +  M) 


! 

I 


4  - 

=  24  +  32u  +  8u2 

2 

♦  I8v  +  24uv  +  6u  v 

+  3v2  +  4uv2  +  v2v2,  (4) 

The  resulting  brightness  values  of  the  image  are: 

P-» 


To  ensure  non-negativity  of  the  image  the  roots 

a  b„  must  be  positive  reals,  or,  if  complex,  occur  in 
r  ,  a 

complex  conjugate  pairs  with  positive  real  components. 

Fried  has  shown  that  the  autocorrelation  function 
of  an  image  which  is  completely  factorable  is  not 
unique,  i.e.,  there  are  many  non-negative  images,  all 
completely  factorable,  which  produce  the  same  auto¬ 
correlation  function.  These  images  are  found  by 
replacing  any  combination  of  the  roots  by  their 
reciprocals,  a"1’  b”1.  For  example,  for  the  factors  (u 
+  a1)(v  +  b^)  the  complete  set  of  images  would  be  gen¬ 
erated  by 

(u  +  a  ^ )  ( v  +  b.j) 

(u  +  a”1 )  (v  +  b^ 

(u  +  a^) (v  +  b”1 ) 

(u  +  a"1 ) (v  +  b"1 )  . 


(5) 


r 
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For  a  polynomial  of  R  u  factors  and  S  v  factors  an 
image  of  dimensions  (R  +  1)  by  (  S  +  1)  results.  By 
flipping  the  roots  2^+T  images  can  be  generated,  all 
having  the  same  autocorrelation  function.  One  of  these 
is  the  original,  the  other  its  1800  rotation  so  that 
2s+t  -  2  =  2S+^"1  ambiguous  images  can  be  formed. 

It  should  be  noted  that  the  image  developed  from 
this  method  are  very  constrained  in  terms  of  degrees  of 
freedom.  The  images  are  really  the  product  of  two  one 
dimensional  images,  each  extended  in  a  direction 
orthogonal  to  its  variable.  Let  f(p)  be  the  one  dimen¬ 
sional  image  obtained  from  the  expansion  of  the  u  roots 
and  f(q)  the  one  dimensional  image  obtained  from  the 
expansion  of  the  v  roots.  the  two  dimensional  image 
f(p,q)  is  then 


f ( p ,q)  =  f ( p)  f ( q )  . 


(6) 
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PRELIMINARY  IMAGES 

A  program  was  written  by  B.  Fahy  of  the  visibility 
laboratory  which  allows  specification  of  the  u  and  v 

factors  or  roots,  which  roots  are  to  be  flipped,  and 

computes  the  resulting  image,  normalized  so  that  the 
maximum  value  is  10,000.  Results  are  shown  in  Tables  1 
to  25  for  a  variety  of  root  distributions.  In  all 
cases  the  image  brightness  values  decrease  monotor.i- 
cally  from  the  maximum  value.  As  the  roots  increase  ir: 
value  from  unity  the  images  become  smaller,  that  is, 

fall  off  more  rapidly.  They  also  shift  towards  an  edge 
and  become  more  asymmetrical.  Flipping  half  of  the 

roots  changes  the  location  and  the  symmetry.  In  gen¬ 
eral  the  images  produced  are  uninteresting. 
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FACTORS  AS  POINT  SPREAD  FUNCTIONS 

The  process  of  multiplying  a  polynomial,  whose 
coefficients  represent  an  image,  by  a  factor  (u  +  ar) 
is  equivalent  to  convolving  that  image  with  a  pair  of 
delta  functions.  For  example  the  polynomial 

=  A0  +  A-jU  +  A2u^  +  (7) 

when  multiplied  by  (  u  +  ap  )  will  consist  of  two  com¬ 
ponents.  The  first  will  be  the  original  polynomial 
with  each  coefficient  multiplied  by  ar  .  The  second 
component  will  be  the  original  polynomial  with  the 
power  of  each  variable  increased  by  one,  which 
corresponds  to  a  right  shift  of  the  image  associated 
with  the  polynomial  coefficients.  In  terms  of  the 
image  we  have  convolved  the  original  image  with  pair  of 
delta  functions,  one  at  the  origin  with  a  weighting  of 

an  and  the  other  shifted  to  the  right  by  one  pixel, 
with  a  weighting  of  unity. 

If  the  polynomial  is  multiplied  by  ( +  aj  where 
6  is  an  integer,  the  effect  on  the  one-dimensional 
image  associated  with  the  polynomial  is  to  convolve  it 
with  a  pair  of  delta  functions,  now  separated  by  6  pix¬ 
els.  This  allows  more  interesting  images  to  be  con¬ 
structed,  that  is,  images  with  more  structure.  In 
addition,  the  ambiguous  images  can  be  significantly 
different  from  one  another.  This  is  illustrated  in 


Fig.1  which  is  a  graph  of  two  one-dimensional  funci- 
tons.  The  change  in  going  from  one  image  to  the  other 
was  caused  by  the  flipping  of  one  root. 


The  convolution  concept  leads  to  a  less  con¬ 
strained  method  of  generating  ambiguous  images.  It 
will  now  be  shown  that  any  arbitrary  one-dimersional 
function  can  be  used  to  produce  a  pair  of  ambiguous 
images  by  convolving  the  image  with  a  pair  of  delta 
functions.  The  initial  function  does  not  need  to  be  a 
factorable  function. 

Let  f(p)  be  a  general  or.e  dimensional  real  func¬ 
tion.  Convolving  f(p)  with  the  delta  function  pair 
produces  a  new  function  g(p),  which  is  the  original 
function  weighted  by  (a)  plus  the  original  function 
shifted  by  6. 

g(p)  =  a  f(p)  +  f(  p-6)  (8) 

The  autocorrelation  function  of  g(p),  C(i),  is 

2g(p)  g(p-i) 

C(i)  =  £ - 5 -  (9) 

Zg( p) 2 

We  substitute  Eq.(3)  into  Eq .  ( 9 )  and  use  the  relation¬ 
ship  of 

3f(p)f(p-i)  =  2f(p-6)f(p-6-i)  OO) 

P 


to  obtain 
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C ( i)  = 

( 1+a2)Sf( p) f( p-i)  +  aJ[f(p-6)f(p-i)  ♦  f( p) f( p-6-i) ] 

- P _  (11) 

(1+a2)  2f2(p)  +  2a2f ( p) f( p-6) 

If  we  multiply  top  and  bottom  of  Eq (11)  by  1/a2,  which 
leaves  C(i)  unchanged,  the  resulting  equation  is  ident¬ 
ical  to  that  obtained  by  replacing  a  by  1/a.  In  other 
words  convolving  the  function  f(p)  with  the  two  point 
spread  functions 

a 

1 

P 

0  6 


1 

1 

II _ 

- y 

0  6 


produces  images  with  the  same  autocorrelation  function. 
Again  f (p )  can  be  any  function,  it  need  not  be  derived 
from  a  factorable  polynomial. 
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Figure  2  shows  a  more  complex  set  of  ambiguous 
images.  The  starting  image  was  an  array  of  seven  dots, 
arranged  in  the  pattern  found  on  a  domino.  The  factor, 
_i.e.,  the  horizontal  shifts  and  weightings  used  to  gen¬ 
erate  the  ambiguous  images  are  shown  on  the  figure. 
The  four  images  have  identical  autocorrelation  func¬ 
tions.  This  was  verified  computationally.  The  four 
images  have  the  same  basic  structure  in  terms  of 
numbers  and  locations  of  dots  but  the  relative  bright¬ 
ness  of  the  dots  are  quite  different.  The  convolving 
process  could  be  done  in  the  vertical  direction  to  pro¬ 
duce  even  more  complex  ambiguous  images. 


11 


DISCUSSION  AND  CONCLUSIONS 

The  work  shown  here  has  demonstrated  that: 
a)  Factorable  polynomials  using  unity  power  factors, 
such  as  (u  +  an),  result  in  smooth  images  which  change 
little  apart  from  spatial  shifts  when  ambiguous  images 
are  generated;  b)  When  greater  than  unity  power  factors 

y/nx 

are  used,  such  as  (u/v  '  +  an),  structured  images  can 
result  and  the  ambiguous  images  can  be  quite  different 
from  one  another;  c)  The  concept  of  multiplication  of  a 
polynomial  by  a  factor  is  equivalent  to  convolution  of 
the  image  represented  by  the  coefficients  of  the  poly¬ 
nomial,  with  a  pair  of  delta  functions.  This  concept 
leads  to  the  result  that  the  starting  function  for  gen¬ 
erating  ambiguous  images  need  not  be  a  factorable  poly¬ 
nomial  but  can  be  any  non  negative  function.  That  is, 
convolution  of  any  starting  function  with  appropriate 
pairs  of  delta  functions  can  lead  to  different  images 
which  have  the  same  autocorrelation  function.  An  exam¬ 
ple  of  four  ambiguous  images  using  the  technique,  shown 
in  Fig.  2,  shows  that  while  basic  shape  and  structure 
was  the  same  for  all  images,  considerable  variation  in 
the  internal  brightness  values  occurred  between  the 
images . 

Images  which  are  ambiguous  are  very  specialized, 
at  least  as  generated  by  this  technique.  Such  images 
may  not  occur  frequently  in  practice.  However,  if 
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image  reconstruction  is  made  from  the  autocorrelation  function, 
the  possibility  of  the  reconstructed  image  being  ambiguous  should 
be  considered. 

At  the  present  there  apparently  is  no  theory  available  which 
allows  determination  of  whether  a  given  autocorrelation  function 
is  unambiguous,  that  is,  whether  the  autocorrelation  function 
could  have  been  produced  by  only  one  non  negative  image. 
Instead,  the  theory  deals  with  the  reverse  of  the  problem: 
methods  of  generating  images  which  have  either  ambiguous  or 
unambiguous  autocorrelation  functions. 

Perhaps  the  approach  for  now  is  to  rely  on  the 
reconstruction  techniques  themselves  for  answer  to  the  question 
of  uniqueness.  For  a  wide  class  of  images,  with  and  without 
noise,  do  the  techniques  provide  reconstruction  of  the  image 
within  limits  set  by  resolution  of  the  optical  system  and  sensor 
noise?  Comparisons  should  be  made  to  other  image  processing 
approaches . 
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Figure  2.  Ambiguous  images  having  the  same  autocorrelation 
function.  Starting  function  was  a  pattern  of  seven  dots  arranged 
as  on  a  domino.  Let  f(u,v)  be  the  polynomial  (not  factorable) 
representing  the  starting  dot  pattern.  The  images  are: 


Top  Left 
Top  Right 
Bottom  Left 
Bottom  Right 


f(u,v)(u7+2)(u14+2) 
f(u,v)(u7+2)(u14+.5) 
f(u,v)(u7+.5)(u14+2) 
f(u,v)(u7+.5)(u14+.5) . 


o  o 
o  o 


-H  h-  14  pixels 


Original  pattern: 


TABLES  1  to  25 


These  tables  represent  the  brightness 
values  of  two  dimensional  images  made 
from  completely  factorable  polynomials 
whose  factors  are  of  unity  power. 


t 


i 


TABLE  2 


TABLE  3. 

1/2  of  roots  s  2. 
1/2  of  roots  =  1/2. 


TABLE  4. 


END  IMAGE 


1/2  of  roots  = 


end  image 
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TABLE  8 


TABLE  9 


TABLE  10. 


TABLE  11 


END  IMAGE 


TABLE  12. 
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TABLE  13. 


Linearly  Increasing  roots,  1  to  10. 
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Linearly  Increasing  roots,  1  to  10, 

«t  half  fllooed.  Compare  to  TABLE  14. 


TABLE  16. 

Linearly  Increasing  roots,  1  to  10, 


TABLE  17. 


Linearly  increasing  roots,  10  to  100, 
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TABLE  19. 

Linearly  increasing  roots,  10  to  100, 
icond  half  flipped.  Compare  to  TABLE  17. 
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Exponentially  increasing  roots,  /?  steps. 
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END  IMAGE 


TABLE  21 . 

Exponentially  increasing  roots,  /?  steps. 
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Exponentially  increasing  roots,  /2  steps. 


o  o  o  ~  I  ro  n  o  O'  I  o  n  o 


O  O  "4 

o» 


O  O  N  0* 

O'  I  •€  ID 


O  O  N  <»■ 

-  ru 
ru 


o  o  o  T 
oi  o 
01 


O'  01  <0  O  O' 

■o  n  co  o  n 

'T  -<  -i  K 

co  n  n 


O  O'  01  N  O 

O  <0  O  01  N 

O  ♦  S  CO 

o  co  n 


000f0-»0»00l0*«0) 

-'diviv^nmo^ 


o  o  n  ^ 

♦ 


N  N 
'O  O' 
N  + 


OOON'OO*** 

•  ♦con 

01  01 


oi  oooonoioi'' 

o'  ~  01  ~ 

< 


<c  o  o  o  o 

LU 

a: 


»-  o  o  o  o 

D 

O 

GO 

< 


L 


jo  o< 

30  Q( 


END  INAGE 


r 


t/J 

CO 

Q. 

CM 

Q) 

4-> 

LU 

</1 

_t 

CQ 

04 

< 

X 

l-~ 

*» 

O 

+J 

44 

O 

Ci 

o 

L. 

• 

L. 

«J 

Q. 

CM 

a» 

6 

c 

o 

ui 

*r- 

<_> 

i/t 

CO 

10 

<c 

01 

• 

H- 

u 

u 

•S 

c 

Cl 

a. 

>> 

r» 

4- 

<o 

4. 

•r* 

+4 

«0 

c 

a/ 

c 

+4 

o 

VI 

a 

s- 

X 

UJ 

u_ 

t 


V 


-  39 


* 

ui 

CP 

< 

E 


Q 

< 

UJ 

a 

a 

i- 

\r 

O 

09 

< 


o  o 


o  o 


o  o 


o  o 


o  * 


ft»  n 
oo 


«r  o 
m 


o i  ni 
o> 


O  CVI 

m 


O  (Vi 


o  o 


o  o 


o  o 


o  o 


o  n 


4r  (*■ 


N  O 
▼  O 
CD 

n 


*  (VI 

«-<  ~i 

<n  -* 
-o 


to  <0 
cvi  n 
0“  n 
m 


%  8 
ni  o* 


(VI  (Vi 
(Vi  0* 


o  rvi 


o  o 


o  o 


o  o 


«  n 


(vi  pj 

H  S 

n  ~ 


2  £ 
«0  O' 

fi  o 
m  n 


o  n 
o  o 
o  in 
o  tr- 


cvi  o 
^  o 

-i  *r 
*o  n 


*  cvi 

to  00 


o  n 
to  oo 


*  (VI 


o  o 


o  o 


o  o 


o  o 


o*  o 


n  n 

N 


cvi  to 

H 

n 


«  n 

O' 


fs  o 

* 


•r  o 


o  o 


o  o 


o  o 


o  o 


o  o 


o  o  o 


o  o 


o  o 


o  o 


o  o 


o  o 


o  o 


o 

o 

ni 

^4 

n 

# 

o 

o 

•o 

m 

ni 

* 

O 

o 

03 

ni 


o 

o 

to 

ni 


o 

o 

o 

n 


o 

o 

is 

r 


■s 

s 

1 

8 

8 


1 

» 

* 

1 

1 

O 

o 

1 

1 

o 

o 

1 

o 

o 

1 

o 

o 

1 

1 

<t 

i 

1 

o 

<0 

1 

o 

o 

H 

o 

o 

o 

o 

o 

o 

ni 

..  ni 

n 

n 

m 

n 

ru 

ni 

•a 

<0 

o 

o 

6 

6 

i 

UJ 

i 

C9 

o 

o 

< 

n 

m 

E 

ru 

eu 

H 

T+ 

Jl 

■  1. 

a 

6 

o 

z 

UJ 

no  ui( 
»>o  tS 

<18  8 

!U  5 


Exponentially  Increasing  roots*  x2  steps. 
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